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The Millisecond Radio Sky: Transients from a Blind Single 
Pulse Search 
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ABSTRACT 

We present the results of a search for transient radio bursts of between 0.125 and 32 
millisecond duration in two archival pulsar surveys of intermediate galactic latitudes 
with the Parkes multibeam receiver. Fourteen new neutron stars have been discovered, 
seven of which belong to the recently identified "rotating radio transients" (RRATs) 
class. Here we describe our search methodology, and discuss the new detections in 
terms of how the RRAT population relates to the general population of pulsars. The 
new detections indicate (1) that the galactic z-distribution of RRATs in the surveys 
closely resembles the distribution of pulsars, with objects up to 0.86 kpc from the 
galactic plane; (2) where measurable, the RRAT pulse widths are similar to that of 
individual pulses from pulsars of similar period, implying a similar beaming fraction; 
and (3) our new detections span a variety of nulling fractions, and thus we postulate 
that the RRATs may simply be nulling pulsars that are only "on" for less than a 
pulse period. Finally, the newly discovered object PSR J0941-39 may represent a link 
between pulsars and RRATs. This bizarre object was discovered as an RRAT, but in 
follow-up observations often appeared as a bright (~10 mJy) pulsar with a low nulling 
fraction. It is obvious therefore that a neutron star can oscillate between being an 
RRAT and a pulsar. Crucially, the sites of the RRAT pulses are coincident with the 
pulsar's emission, implying that the two emission mechanisms are linked, and that 
RRATs are not just pulsars observed from different orientations. 

Key words: stars:pulsars — stars:neutron stars — methods: data analysis 



X 



1 INTRODUCTION 

The transient radio sky has revealed many classes of ex- 
otic objects since the discovery of pulsars in late 1967 
l|Hewish et al.|[l968l ). While the first ~50 pulsars were dis- 
covered by searching for their single-pulse emission, the ad- 
vent of filterbanks and digital signal processing searching in 
the sub-second radio sky to use more efficient and success- 
ful allowed the use of Fourier search techniques, which are 
more sensitive to periodic objects than single pulse searches 
and have since permitted the discovery of more than 1500 
pulsars. However, the recent disco very of sparsely burstin g 
rotating radio transients (RRATs, iMcLaughlin et al.| [2006). 
and a one-off bur st tentatively identi fied as having an ex- 
tragalactic origin l|Lorimer et al.ll2007h has spurred renewed 
interest in the search for single-pulse emission after 30 years 
of Fourier transform searches dominating surveys for pul- 
sars. 

The more recent history of single-burst searches corn- 
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menced with the search of I Amy et al.l (|l989l ) for impulsive 
signals between 0.001 and 800 ms using a transient event 
monitoring system at 843 MHz at the Molonglo cylindri- 
cal telescope. While they found only pulsars and ambiguous 
events, they were able to put an upper limit on the event 
rate of 1 — 25 ms duration events of 0.017 events s -1 sr -1 
for pu lses with energy density > 10~ 28 J m~ 2 Hz -1 . iNice] 
(1999) also recognised that some pulsars might be unde- 
tectable in periodicity searches and subsequently re-searched 
data from a 68 deg 2 area of sky along the galactic plane 
for single pulses with duration between 0.5 — 8 ms, reveal- 
ing one new, slow (2.3 s) pulsar. The survey which discov- 
ered RRATs (|McLaughlin et al.l [20061 , hereafter M06) was 
the first large-scale examination of a significant data set for 
impulsive bursts since the early pulsa r surveys, while the 
search of the Magellanic Cloud surveys l|Lorimer et al.ll2007t) 
and a tran sient search with Arecibo telescope ( Deneva et al.l 
lsubmitted[ hereafter D0 9) and a deeper search of the M06 
data ( Keane et al.ll2009l . hereafter K09) followed using sim- 
ilar techniques. 

There are several known sources of fast transient ra- 
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dio emission which we expect to detect as single pulses on 
sub-second time scales. Periodic rotators which have irregu- 
lar emission properties include the aforementioned RRATs, 
which are a form of rotating neutron star, emitting pulses 
very sparsely at several to many times the rotation pe- 
riod of the star. At present it is unclear from the few ob- 
served RRATs how they relate to the neutron star popu- 
lation as a whole, or how they may relate to other sub- 
populations of neutron stars such as anomalous X-ray pul- 
sars or magnetars. It is possible that they represent only a 
phase through w hich all pulsars undergo when they reach a 
suffic ient age (e.g. lMcLaughlin et al . 2009; Keane fc Kramer! 
2008), or that are simply distant, low-l uminosity pulsars 
with e xtreme pulse amplitude modulation (jWeltevrede et ahl 
l2006bl ). Such distinctions are important for reconciling the 
so-called "birthrate problem"; it has been recently noted 
that the net neutron star birthrate calculated from appar- 
ently disparate neutron star populations (X-ray dim iso- 
lated neutron stars (XDINS), magnetars, RRATs, pulsars) 
exceeds the rate of the core-collapse supernovae that are be- 
lieved to produce them (|Popov et a l. 2006; Keane fc Kramerl 
2008). Recent works are endeavouring to perform analyses 
that make the distinct ion clearer from pulse timing (e.g. 
iMcLaughlin et al.ll2009l ). 

High pulse-to-pulse modulation and periodic nulling 
represent types of signal variance observed in the single pulse 
emission of some pulsars (e.g. iHerfindal fc Rankinl 120071 ; 
iKramer et al.|[2006l ). Fourier transform searches become less 
sensitive to such sources due to the objects' erratic emission 
and the often short duty cycles and long rotation periods of 
the stars. Such neutron star-related phenomena promise in- 
sights into the localisation and power source of radio pulsar 
emission. Single pulse searches have the ability to recog- 
nise these sources, and may serve to reveal a bridging pop- 
ulation or define a distinction between regular pulsars that 
emit steadily and the seemingly more rare "part-time" pul- 
sars that only emit detectable emission some fraction of the 
time, revealing whether RRATs are (in physical mechanism) 
just extreme examples of the nulling and modulation seen 
in some pulsars. 

Theoretical arguments predict fast transients from a 
variety of non-neutron sta r origins, such as coalesci ng sys- 
tems of relativistic o bjects dHansen fc Lvutikovll200ll ) , evap- 
orating black holes llReesl Il977l ). or supernova events (e.g. 
IColgate fc Noerdlinger]|l97ll ) . The aforementioned detection 
of a single, apparently ext r agalac tic, 5 ms and 30 Jy burst 
detected bv lLorimer et all (|2007h indicated that the detec- 
tion of such sources in other galaxies might be possible. 
While the Lorimer et al. detection was the primary motiva- 
tor for performing our search, fourteen other discoveries as- 
sociated with intermittent neutron star phenomena resulted, 
and these are presented and analysed below. No further de- 
tections were made of bursts attributable to an extragalactic 
origin. Fifteen seemingly non-astrophysical, frequency- swept 
bursts were discovered with a similar frequency-dependent 
sweep rate (or dispersion measure) to the Lorimer et al. de- 
tection (except that they occurred in all 13 instead of just 3 
beams of the Parkes multibeam receiver), however a deeper 
analysis and the interpretation of those signals will be pre- 
sented elsewhere (see Burke-Spolaor et al., in prep). 

This paper presents the detections from a blind search 
for single pulses in more than 900 hours of archival data from 



the pu lsar surveys of lEdwards et all (|200ll ) and lJacobv et all 
(2009). The data is described in section [2] and the search 
methodology in section [3] The fourteen new objects and 
other results uncovered by this search are presented in sec- 
tion [4] The implications of these detections are discussed in 
section [S] 



2 DATA SAMPLE 

The data used in this search were the southern-sky 
(-100° < I < 50°) intermediate (5° < |6| < 15°) and high 
(15° < |6| < 30°) galactic latitu de surveys of lEdwards et all 
j200ll ) and IJacobv et al l (|2009l ). respectively; here b and I 
represent galactic latitude and longitude, respectively. Both 
surveys employed the same analogue 1-bit filterbank as a 
back-end to the 20cm multibeam receiver system on the 
Parkes Radio Telescope. The multibeam receiving system 
allows the simultaneous observation of 13 positions on the 
sky w ith approximately 30 arcmi nutes between beam posi- 
tions (|Stavelev-Smith et al.|[l996T ). The twelve outer beams 
are arranged in two concentric sets of six (see the upper left 
panel of Fig. [1} , with a full- width half power beamwidth of 
each beam at 14 arcminutes. Due to the sharp sensitivity 
fall-off within 24 arcminutes of each beam's pointing centre, 
a pointlike astrophysical signal will generally be detected in 
no more than four beams of the receiver, and unless partic- 
ularly bright, usually just one or two. 

Each survey sampled 96 channels across a 288 MHz 
band centred on a frequency of 1372.5 MHz. The data were 
recorded after the polarisations were summed and then sam- 
pled at 125 fis. Each 13-beam pointing of the telescope lasted 
4.4 minutes. The data total 12,090 beam hours (there are 
13 beams), taken over the years 1998-1999 (Edwards) and 
2001-2002 (Jacoby). The total number of pulsars detected 
in the periodicity searches of these surveys were 170 for Ed- 
wards et al. and 62 for Jacoby et al, giving a total of 230 
independently detected pulsars. 

All follow-up data presented here were taken with the 
new 2-bit digital filterbank system currently installed at 
Parkes Radio Telescope that use CASPER tools and is 
known as BPSR. Data with this system were taken with 
1024 channels across 400 MHz centred at 1381.8 MHz, and 
sampled at 64 /is. Filters blocking emission from a geosyn- 
chronous satellite visible by the telescope render frequencies 
above 1520 MHz unusable. 



3 SINGLE BURST SEARCH PROCESS 

The most significant difficulty in single pulse searches is al- 
laying the millions of spurious candidates which arise from 
local man-made and natural radio frequency interference, 
such as lightning. Our search employs various stages of auto- 
matic radio frequency interference (RFI) elimination before 
candida tes are manually assess e d. We generally followed the 
steps of lCordes fc McLaughlin! (|2003l ). however with adjust- 
ments in the techniques used for RFI rejection, matched fil- 
tering, and visual inspection. The main discriminant used in 
single pulse searches is the fact that astrophysical transients 
exhibit a well-defined frequency-dependent delay caused by 
the propagation through multiple ionised plasma screens in 
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SNR = 17.74 RA = 1S2553.54 

DM = 40.78 Dec = -332217.39 

tscr = 32.00 glat = not provided 




Figure 1. A detection plot for a burst from J1825— 33. Upper left: Time series from each of the 13 beams dedispersed at the DM of 
brightest detection; each field plots the power versus time, normalised by the RMS of the time series. All fields plot the same abscissa and 
ordinate range, and the beam (and sample in the time series) of brightest detection is indicated by a red star. In the other beams, the 
sample of peak detection is marked by a vertical dotted line, and the detected width of the pulse is marked by a horizontal bar. Lower, 
left: Time series for the maximum-SNR beam dedispersed at the DM of brightest detection, stacked with the time series from one DM 
trial higher ("lup"), one trial lower ("ldn"), and the DM = time series. As in the upper left panel, the ordinate gives normalised power 
(SNR), while the abscissa gives time since the beginning of the file (in seconds). Bottom RIGHT: A false colour image of the received 
power as a function of frequency (in MHz) versus time since the beginning of the file (in seconds), revealing the frequency-dependent 
structure of the detected pulse in the beam of brightest detection. Upper. RIGHT: Signal-to-noise ratio as a function of DM trial and 
boxcar filter size, and best-fit information about the detection (DM, SNR, pointing position of the beam of brightest detection, and 
boxcar width "tscr"). Pulses with resolved complex structure can show multiple peaks and other structure in th eir SNR vs. DM curve, 
thoug h genuinely dispersed pulses (and the curve above at the correct DM±5) generally obey the function given bv lCordes &: M cLaughlin 
1)20031) . 



the interstellar medium, visible as a time lag across the 
bandwidth of our instrument. The dispersion measure, DM, 
measures the integrated line-of-sight column density of free 
electrons in units of pccm -3 . Locally generated signals gen- 
erally do not appear dispersed (DM ~ 0), while genuine 
bursts have a decrease in measured significance as the trial 
dispersion increases from the true dispersion of the source. 
Periodic sources will show multiple detections at the same 
dispersion measure. The radio interference "forest" local to 
any given astronomical observing site gives rise to the need 
for a set of thresholds that remove broad, zero-dispersion sig- 
nals. We therefore have increased confidence in candidates 
at higher dispersion. 

Our detection thresholds were estimated by manual in- 
spection of the data and trial runs on a sub-sample of data to 
characterise the typical local RFI environment of the Parkes 
Telescope. These tests indicated that events separated by at 
least 3.75 ms (30 data bins at the native time resolution) 
could be considered independent, that signals peaking in 



signal-to-noise ratio (SNR) below DM = 5 were most likely 
RFI or statistically insignficant, and that setting a search 
threshold at six times the RMS noise of the dedispersed 
time series was sufficient to filter out the majority of can- 
didates which were due to statistical fluctuations or narrow 
band interference. 

The search process followed several steps. At each trial 
DM, after dedispersion the frequency channels were summed 
and a search for SNR^ 6 pulses was performed of the com- 
plete 4.4 minute observation. An increased sensitivity search 
for broader pulses was then performed, using a progressive 
box-car smoothing function ranging 2 1 ,2 2 ,...,2 8 time sam- 
ples, corresponding to a maximum filter width of 32 ms. 
The 32 ms maximum width was chosen to significantly de- 
crease the ratio of local to astrophysical candidates (as de- 
termined by tests on a data subset), and to preserve our 
sensitivity to impulses of low dispersion measure (the total 
delay in milliseconds across our band is roughly equivolent 
in value to the dispersion measure of an impulse in units of 
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pccm -3 ). We have decreased sensitivity to pulses of half- 
maximum width, w, outside the range 0.125 < w < 32 ms. 
Additionally, the one-bit nature of the data made it diffi- 
cult to discover pulses with width less than a sample pe- 
riod, as we only had ^ 96 functioning one-bit channels, and 
an RMS that meant it was difficult to have much confi- 
dence before summing samples. We searched in DM from 
< DM < 600 pccm -3 , which for the galactic latitudes of 
the survey is sufficient to probe out to at least 50 kpc except 
in the regions a few tenths of a degree from I — 0, |fc| — 5. 

For each beam, a first-pass RF1 rejection was performed 
by doing temporal coincidence matching for all trial disper- 
sions and matched filter steps; in this way we were able to 
reconstruct multi-plane information about singly-detected 
events. A further temporal coincidence search of all remain- 
ing events was done in the 13 beams, conservatively reject- 
ing all candidates which had detections coincident in more 
than nine. After this step, we have a full event archive in 
which are encoded the widest duration statistics for detected 
beams, and information about the beam, DM, and smooth- 
ing filter of strongest detection. Event archives which con- 
tained less than three events, and those which had a peak in 
signal-to-noise at less than DM = 5, were rejected. This step 
significantly reduced the number of high dispersion candi- 
dates which were due purely to narrow-band RFI and small 
number statistics. The information in the remaining event 
archives were used to make the plots used in manual candi- 
date assessment. 

The manual candidate ranking involved the viewing of 
each candidate using the tool displayed in Figure [T] which 
shows a single pulse detection from one of our new dis- 
coveries, and contains the information with which one can 
quickly discriminate a true burst from a false detection. 
During candidate viewing, we simultaneously queried the 
ATNF pulsar catalogue for pulsars within two degrees of 
the brightest beam, such that we were able to easily deter- 
mine DM and position incidence, and attribute detections 
to previously known pulsars. Non-RFI detections not asso- 
ciated with known pulsars were considered new transient 
candidates. 

Assuming that detections which show repeated pulsa- 
tions are part of the neutron star population, we fit a rota- 
tion period to each object which has emitted at least three 
detectable pulses. Pulse periods were derived by first fitting 
periodicities by eye to the detected pulses. The data were 
then broken into sub-integrations at the estimated period, 
and sub-integrations containing on-pulses were extracted 
and used to do an optimised automatic search in period 
and DM space using the pdmp program distributed with the 
freely-available PSRchive software suit^E Where follow-up 
data was available and the observations were of higher SNR, 
those data were used to compute the fits. 



4 RESULTS 

Our blind search yielded the detection of several thousand 
single pulses with which we could associate 91 pulsars, that 
were already previously detected through Fourier transform 
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search techniques by Edwards et al. (2001) and Jacoby 
(2005), indicating that roughly 40% of the pulsars detected 
in their periodicity searches have one or more detectable 
pulses per 4.4 minute pointing. The pulsars we detected do 
not appear to have any selective distribution in age or mag- 
netic field strength, have periods ranging 0.089 to 8.5 sec- 
onds, and appear to have primarily populated the "island" 
of standard pulsars in the period-period derivative diagram 
with periods between ~0.1 and 3 seconds. 

The new objects detected in this search include two 
sources which emitted only one (albeit very convincing) 
pulse over the time we have observed them, and 12 sources 
from which multiple pulses were observed. All of these new 
sources are interpreted to be erratically emitting neutron 
stars. The best-fit parameterisations for each source are pre- 
sented in Table [TJ and a more detailed description for each 
is given below. The implications of the new detections are 
discussed in 



4.1 New Transient Sources 

The properties of the pulsed sources presented here give 
no indication that they are associated with a non-neutron 
star origin. Below we describe the properties unique to 
each new object, and note the amount of follow-up time 
which has been dedicated to each. Where emission was 
detectable in a Fourier search, the intermittency ratio 
(r = SNRsin g ic P uisc/(SNR) Fouricr , as given by D09) is pro- 
vided. 



4.1.1 PSR J0735-62 

This source emitted several pulses with SNR > 7 and many 
pulses below a SNR of 6 in the discovery pointing, and af- 
ter a period was determined from the brightest bursts, the 
source was visible when the data were folded at the manually 
measured period with a SNR of just below six. The source 
appears to experience severe scintillation of its signal due 
to its low dispersion measure, and it is likely that because 
of this the six-minute tracks on this source over 93 minutes 
of follow-up observations showed no detectable emission in 
single pulses, however showed a weakly visible signal when 
the data were folded at the nominal period and DM. Three- 
period averaged pulses from this pulsar are shown in Figure 
|3(a)| Despite its steady emission, the long period and narrow 
duty cycle prevent it from being d etected in standard Fourie r 
searches of short duration (e.g. iDeneva et al.l Isubmittedl) . 
The intermittency ratio of this object is r = 1.45. 



4.1.2 PSR J 0923-31 

This source emitted only one pulse in the discovery observa- 
tion, which appeared at nine times above the system noise 
and with complex profile structure with a bright, narrow 



main component. This pulse is shown in Figure 4(a) and 



appears to have a faint, time-resolved, double peak. Sub- 
sequent pointings, totalling 30.9 minutes, revealed no fur- 
ther emission from this location. Further pulse detections 
are required to rule out local interference as the source of 
the burst, however the appearance in only one beam of the 
multibeam receiver and frequency sweep that closely follows 
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Figure 2. A galactic map of the pulsars originally detected in the data, plotted with the new detections presented in this work. 
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Table 1. Data for the 14 new sources detected by this search. Columns: (1) Name based on J2000 coordinates (# signifies only one burst 
detected, while f denotes an object that was Fourier-detectable in the original pointings or in follow-up data; * see t|4. 1.31 for special 
remarks on J0941-39); (2,3) galactic longitude (I) and latitude (6); these coordinates have an error of approximately ±15 arcminutes due 
to the poor positional localization of the multibeam rec eiver; (4) best-fit disp ersion measure and error; (5) distance to the object based 
on the DM and NE2001 galaxy electron density model l|Cordes fc L azio 2002); (6) the pulsation rate based on all data taken on object 
position to date, where the pulsar was not detected above a SNR of 5 in the folded data. An asterisk indicates an object for which only 
the 4.4-minute detection pointings have been used to compute this number; (7) the best fit period to all detected pulses with error on 
the last digit in parentheses; (8) integrated pulse width measurement at 50% of maximum after the stacking of all detected single pulses. 
For PSR J 0923— 31, we quote the narrow component width and the total pulse width; (9) pulse width measurement at 50% of the peak 
of the brightest detected pulse. PSR J0923-31 specified as in column 8; (10) peak flux density of brightest detected single pulse. 



the cold plasma dispersion relation gives a strong implica- 
tion that the source is astrophysical. Based on a 10% duty 
cycle and a pulse width of 213 ms, if this object is a pulsar 
its estimated period should be > 2 seconds. 

4.1.3 PSR J0941-39 

This pulsar was discovered after the detection of five sin- 
gle pulses to which there was no clear periodicity. Follow- 
up observations of the source (totalling more than 3 hours 
spread over approximately one-half year) later revealed ex- 
traordinary behaviour, sometimes appearing with sporadic 
RRAT-like emission at a rate of ~90-110 per hour, and at 
other times emitting as a bright pulsar with strong pulse- 
to-pulse modulation and quasi-periodic nulling (see Fig. 



O. The pulsar has a triple-peaked profile and broad duty 
cycle, with complex sub-pulse drift, mode-changing, and 
longitude-dependent modulation when in an "on" state (it is 
these behaviours which made single-pulse periodicity fitting 
difficult). A timing program is underway on this pulsar at 
Parkes Telescope, and observations of the properties of this 
object will be presented in detail in a forthcoming paper. 

4.1.4 PSR J1129-53 

This source was originally detected as two single bursts oc- 
curring at the same DM in independent pointings separated 
by 14 arcminutes. Follow-up observations totalling 15.5 min- 
utes at a position half-way between the original two detec- 
tion positions (it is this position listed in Table [1} revealed 
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Figure 3. Greyscale representations of time (or pulse index, y- 
axis) versus de-dispersed pulse phase — "pulse stacks." These four 
objects were alone in showing Fourier-detectable emission under- 
lying their single bursts. PSR J0735-62 has been folded over three 
rotation periods, PSR J1850+15 over six periods, and the others 
show single pulses. 



10 further pulses which were extracted and used to fit a pe- 
riod and DM. The individual pulses from this object exhibit 
a variety of pulse shapes, with some pulses exhibiting evenly 
double-peaked emission. 



4.1.5 PSR J 1226-32 

Three pulses of SNR > 12 were detected from this object 
during the discovery pointing, and a best-fit to these pulses 
gave us the period and dispersion measure quoted in Table Q] 
When folded at the fitted period, the rotations which did not 
contain pulses of SNR > 12 showed no detectable emission. 
We have not yet had the benefit of follow-up observations for 
this object, however the three dispersed pulses emitted by 
the object show consistency in both dispersion sweep and 
frequency-dependent modulation, indicative of the bursts' 
propagation through the interstellar medium from a com- 
mon origin. 



4.1.6 PSR J1534-46 

This object appears to be intrinsically modulated (and 
nulling) with single pulses appearing at 12-14 times the rota- 
tion period. The modulation power spectrum shows a peak 
at ~ 4 times the pulse period near the centre of the inte- 
grated pulse. The single-pulse emission from this pulsar is 



shown in Figure 3(b) The pulsar was detectable in a Fourier 
search, giving r = 0.66. We have not performed follow-up 
observations of this source. 




i 1 i 1 i 

200 400 600 
Time (ms) 

(a) PSR J0923-31 



800 1000 




i 1 1 1 1 r 
100 200 300 
Time (ms) 

(b) PSR J1610-17 



400 



Figure 4. Greyscale spectra (bottom panel) and dedispersed, 
frequency-integrated timeseries (top panel), showing the single 
pulses of the two objects from which no further emission has yet 
been detected. 



4.1.7 PSR J1610-17 

Only one pulse was detected from this location, however 
there have been no follow-up observations on the object. The 
original detection had a SNR of 7.4, and no other detections 
above SNR = 4 appeared in the 4.4 minute observation. The 



single detected pulse is shown in Figure 4(b) and like PSR 
J0923-31, shows a frequency sweep that agrees well with 
the expected frequency delay due to interstellar dispersion, 
implying a astrophysical, galactic origin of the pulse. 



4.1.8 PSR J 1647-36 

This pulsar appears to exhibit periodic nulling with bursts 
lasting on the order of 1-6 pulse periods spaced by nulls last- 
ing ~ 27 seconds (130 rotations). With each period of "on" 
activity, the bursts arrive at a phase offset by 5 — 10% from 
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Figure 5. Two-minute (205 rotation) pulse stacks exhibiting PSR J0941-39's "RRAT" (off) and "Pulsar" (on) phases f £|4.1.3ll5~3l l. The 
integrated profile shape is shown in the panel the the right. Tracks of the signal-to-noisc ratio of the peak of each pulse with SNR > 3 
are plotted below the pulse stacks. The brightest peaks of the on- and off-phases are comparable, indicating that the "off" state is not 
due to interstellar scintillation. 



the previous burst episode; it is possible that the system is 
either strongly drifting or mode-changing, however follow- 
up observations of this source are necessary to assess these 
hypotheses. 

4.1.9 PSR J1654-23 

The original detection showed 3 narrow (sub-0.5 ms), faint 
(SNR < 7) bursts which allowed the fit reported in Table [1] 
The duty cycle is the smallest observed to date for a pulsar, 
however due to the narrow pulses and weak emission it is 
possible that some pulses were not visible above the noise. 
It is possible that the period reported for this pulsar could 
easily be an integer multiple of the true period, the effect 
of which would be to broaden the duty cycle by the same 
integer value. Follow-up observations of longer duration than 
the discovery pointing and at a faster sampling rate will 
reveal more completely the nature of this pulsar. 

4.1.10 PSR J 1753-38 

The period and DM for this object were derived by extract- 
ing and fitting just the three bright pulses detected in the 
original pointing. Follow-up pointings (totalling 21 minutes) 
revealed further pulses consistent with the flux, period, and 
DM of the original detections. 



strong impulse was present made genuine pulse extraction 
difficult, however three dispersed pulses from the object were 
identified and parameters estimated. We have not yet made 
follow-up observations of this source, however cleaner data 
may reveal more frequent low-level emission. 



4.1.12 PSR J 1825-33 

In the discovery pointing, this pulsar emitted seven bright, 
sequential pulses near the beginning of the observation but 
was not visible thereafter, nor was any emission visible in 
25 minutes of follow-up observations. The profile structure 
of this pulsar has a double peak (Figure Q] shows one pulse 
from this object). 



4.1.13 PSR J1850+15 

This pulsar's discovery pointing contained strong periodic 
interference which probably prevented its detection in the 
original Fourier search (r = 1.04; interference is visible atop 
the fainter pulses of the object in Fig. 3(c) \. The emission 



appears to be steady, with no evidence for clear intermittent 
behaviour, although there is strong pulse-to-pulse amplitude 
modulation. We have not performed follow-up observations 
of this source. 



4.1.11 PSR J 1753-12 

The discovery pointing containing this object was badly af- 
fected by interference. This and the fact that only one very 



4.1.14 PSR J2033+00 

This source appears to be an intermittent pulsar with a high 
nulling fraction and is sometimes detectable in the Fourier 
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domain. There is some evidence of sub-pulse drift; individual 
pulses rarely come from the same phase for more than a few 
rotations (see Figure 3(d) I. The source was detectable in 



single pulses for the entirety of the 25 minute follow-up on 
the source, however became distinctly brighter (and weakly 
detectable in a Fourier search; r = 1.86) in the final 14 
minutes. 



5 DISCUSSION 

5.1 Pulsar transient varieties and the nature of 
RRAT emission 

Regular pulsars are themselves intermittent radio sources, 
showing detectable emission only for the fraction of time 
that their beams sweep Earth. The individual pulses, how- 
ever, are observed to show intrinsic pulse-to-pulse flux and 
shape differences, the origin of which is yet unknown (though 
various theories exist). The variations observed in pulsars 
can be broken down into two main types, being: pulse- 
to-pulse modulation, in which the object emits continu- 
ously, however at an intrinsically modulated amplitude, and 
nulling, during which the pulsar beam is somehow quenched, 
emitting no detectable radiation for a period of time. In 
many pulsars, nulling appears to be a real phenomenon (with 
no underlying emission when in null states), with nulls ob- 
served to occur from zero to 95% of a pulsar's observed 
rotations jWang et al.ll2007t ). Some pulsars appear to null 
for days-weeks, with a p eriod derivative tha t is dependent 
upon the emission state (|Kramer et al.|[2006l ). 

The range of nulling fractions observed in pulsars is eas- 
ily described by a windowing function, in which the number 
of sequential pulses observed at a given time is defined by 
an on-window of some length. In this description, regularly 
emitting pulsars have an indefinite on-window, while pulsars 
with increasing percentages of nulling have on-windows of 
increasingly short timescale, with extreme millers (nulling 
fraction > 75%) becoming more likely to be detected in 
single pulse searches than periodicity searches. Various de- 
grees of pulsar windowing are illustrated in Figure [5] It is 
clear that one could incorporate the RRATs into the gen- 
eral pulsar population in this way, by just extending the 
non-nulling window to a time of length t < P, causing only 
one pulse to be observed at a time. It is a possibility that 
in this model, pulsars start off as continuous emitters, and 
gradually increase their nulling fraction until they become 
RRATs, which ultimately evolve further to emit no radio 
emission. 

The alternative view of RRAT emission, re- 
lated to pulse-t o- pulse modulation, was suggested by 
IWeltevrede et all (|2006bh . who indicated that if PSR 
B0656+14 were observed at several times its distance from 
Earth, the highly variable single pulses (emitting at up to 
420 times the average pulse flux) could emulate the single 
pulse behaviour observed by M06's distant discoveries. 
The pulse-energy-distribution analysis of K09, however, 
makes an argument against this hypothesis (although the 
analysis was not done on "classic RRATs" as we define 
them below). It appears obvious that there will ultimately 
be some fraction of RRATs that belong to this class, but 
deep integrations would be expected to show underlying 



emission; the real issue, therefore, is in what fraction of 
RRATs belong to this class, and if some do not, what are 
the properties of those sources that do not? 

Whether an effect of modulation or nulling, the current 
competing theories of RRAT emission can be broken into 
three possibilities, in which the RRATs represent: 

(1) a distinct class (intrinsic life- long single- pulsing) 

(2) faint pulsars with extreme pulse-to-pulse modulation 

(3) an evolutionary phase of pulsars, e.g. where the mag- 
netic field degrades or the spin period increases, leading to 
increased nulling fractions 

Scenarios (1) and (3) intertwine similarly with questions 
about pulsars observed to null (relating to whether nulling is 
an inherent phase of all pulsars or an effect of the local pul- 
sar environment). The clarification of RRATs' nature not 
only has innate value, but the three scenarios above have 
quite different implications for how the RRAT birthrate in- 
tegrates into that of pulsars, and how they contribute to 
the imbalance of core-coll apse supernova ra t es and neutron 
star birthrates noted bv iKeane fc Kramerl l|2008h . For in- 
stance, in scenario (1) the latent RRAT population may ac- 
tually exceed the pulsar population. Scenarios (2) and (3), 
on the other hand, are in part solutions to the birthrate 
problem, in that the birthrates contributed by RRATs are 
effectively absorbed into those of pulsars. If RRATs are sim- 
ply distant, highly modulated pulsars, the sources detected 
as RRATs are simply abnormalities among the distant, oth- 
erwise undetectable pulsars sitting below the average lumi- 
nosity limit of periodicity surveys; these distant and low- 
luminosit y pulsars are accoun ted for in pulsar birthrate anal- 
yses (e.g. lLorimer et al.lll993h with statistical estimates, us- 
ing luminosity-dependent scaling factors. 

Although the original RRATs of M06 have a number 
of apparent properties which set them apart from regular 
pulsars, as detailed in Sec. 15.31 we define the classic RRAT 
as an object which emits only non-sequential single bursts 
with no otherwise detectable emission at the rotation period. 
When using the term nulling pulsar, we refer to an object 
which emits for some fraction of its rotations, however the 
emission coming in bursts of several sequential pulses. Fi- 
nally, we refer to all objects which emit pulses for less than 
~25% of their rotations as extreme pulsar transients. 

Below we discuss the properties of our new transients 
and compare them with the RRATs of M06 and with the 
properties of the general pulsar population. 



5.2 A breakdown of our discoveries 

What is immediately striking about the 14 objects detected 
in this search is the variety of dormancy timescales, covering 
a broad range of the nulling and modulation phenomena re- 
lated above. For some objects, such as PSR J1850+15, the 
lack of previous detection by Fourier search techniques can 
be explained by the presence in the data of intense periodic 
RFI, which lowers sensitivity to fainter periodic sources. The 
pulsars PSR J0735-62 and PSR J 1850+ 15 appear to sim- 
ply be regular pulsars which were likely missed in previ- 
ous searches due to the modulation and long period of PSR 
J0735-62, and the presence of interference in the discovery 
pointing of PSR J1850+15. 
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Figure 6. Windowing as a descriptor for various emission observed in pulsars, where the average timescale of the window function defines 
the behaviour observed: (a) shows a regular pulsar, (b) is an example of a nulling pulsar, (c) an example of what we term an "extreme 
miller" (e.g. nulling fraction > 75% of pulses are nulled and the on-window typically endures for a few times the rotation period), (d) an 
RRAT, in which the typical on-window has a duration of less than the pulsar period. 



The "classic RRATs" in our search include PSRs J 1129- 
53, J1226-32, J1654-23, J1753-38, and J1753-12. We also 
place the objects from which only one pulse was observed in 
this category; although we have not detected further events 
from the sources, the pulse widths are typical of pulsars and 
the pulsation rates implied by our observations are consis- 
tent with the range of pulsation rates observed in the RRATs 
published by M06, which ranged from one per four minutes 
to one per three hours. We note that the bright single pulses 
emitted by PSRs J1226-32 and J1753-38 are too bright and 
have a pulsation rate too high to represent a long tail of a 
log-normal pulse energy distribution, immediately indicat- 
ing that scenario (2) above (extreme modulation) is not the 
case for all RRATs. This is less clear in the case of PSR 
J1753-12, and we note that our sources may represent a 
mix of highly modulated low-luminosity pulsars, and true 
RRATs. 

Apart from PSR J0941-39, the remaining objects ap- 
pear to be nulling pulsars with high nulling fractions. Only 
two of these objects were in an "on" state for a sufficient 
fraction of observing time to be Fourier-detectable. Addi- 
tionally, it is difficult to say whether the source PSR J1825- 
33 is truly a nulling pulsar, or whether this source had an 
event which briefly activated its outer gap, as we have since 
not detected any further emission from the object. 

PSR J0941-39 does not fit cleanly into any currently 
identified pulsar category, sometimes acting like a single- 
pulse-emitting RRAT, and sometimes like a nulling (though 
not extreme- nulling!) pulsar. This object may play an im- 
portant role in understanding RRAT emission. 



5.3 Exploration of RRATs and the properties of 
our new discoveries 

This search represents an important step in understanding 
extreme pulsar intermittents in several ways. First, we have 
surveyed the largest area to date for such sources in a re- 
gion off the galactic plane, and added the discovery of 12 
highly intermittent pulsars to a small known population of 
such objects. Given that neutron stars are mainly born from 
massive stars in the plane of the galaxy, the high galactic 
latitudes predispose us to sample a population which has 
had a sufficiently high kick velocity and sufficient time to 
reach its current position. Also, the automatic RFI rejec- 
tion and additional candidate visualisation techniques ap- 
pear to have given us a higher sensitivity to low DM single 
pulses, as seen by the fact that M06 discovered no sources 
at DM < 88 pccm~ 3 while we have discovered sources at 
a broad range of DM, with twelve of the fourteen detec- 
tions at less than 88pccm -3 , and with previously known 
pulsars positively identified down to DM = 3pccm~ 3 . Ad- 
ditionally, the latitude ranges of past RRAT searches (K06, 
D09, K09) did not exceed 5 degrees. This increased DM sen- 
sitivity and high, complementary latitude range makes our 
discoveries salient objects for assessing the inherent proper- 
ties of nearby, RRAT-like bursting sources and how neutron 
stars with extreme intermittence relate to the Galactic pop- 
ulation of pulsars. 

Throughout this discussion we analyse all of our sources 
which are sufficiently intermittent to not be detectable in a 
Fourier search — that is, the ten sources not plotted in Figure 
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[3] Where appropriate, we distinguish the analysis of "classic 
RRATs" and that of "nulling pulsars." For reference we also 
include the M06 RRATs and the four D09 objects which, as 
noted by their study, did not show detectable underlying 
regular emission. Additionally, we include the five objects 
which were not noted by K09 to exhibit Fourier-detectable 
emission or clear nulling behaviour. 

The sparse emission and dissimilarity of magnetic 
field stren gth to pulsars which emit giant pulses on fis 
timescales l|Cognard et a.1.1 fl996l ') are what most definitively 
led McLaughlin et al. to the conclusion that their sources 
were a new neutron star population, however several other 
of the observed properties of the original 11 sources also 
distinguish them from the general pulsar population. This 
includes their low duty cycles, high single-pulse luminosities, 
and their tendency to be distant (> 2 kpc, which may have 
partially been due to a selection effect; this particular RRAT 
featur e gives more weight to the claims of IWeltevrede et al.l 
2006b that RRATs may simply be distant, highly modulated 
pulsars). To further our understanding of the properties of 
RRATs, here we explore their various distinguishing features 
in comparison to the extreme intermittents detected in our 
search. 

The galactic distribution of our new objects (of the clas- 
sic RRATs and of our 14 discoveries in general) agree well 
with the distribution of pulsars, concentrated more highly 
towards the Galactic plane and the Galactic centre. Though 
not statistically rigourous with a sample of fourteen objects, 
the galactic z-height distribution (z = d ■ sin(Z)) of our dis- 
coveries above z = 200 pc is well-fit with an exponential 
and a scale height of 300 pc, in good agreement with the 
galactic pulsar distribution. This agreement and the discov- 
ery of RRATs in the high latitude surveys indicate that the 
RRATs are not a young phenomenon confine d to the galac - 
tic plane like anomalous X-ray pulsars (e.g. |Purant| [2006). 
The object of second largest z-height in our sample is the 
RRAT J1610-17 with a scale height of z = 0.86 kpc. With 
a birth at z — pc and a proper motion outward from the 
galactic plane at 100 kms" , the implied age of this object 
is > 10 7 years. 

Figure [7] plots the duty cycles of known pulsars and the 
objects analysed here. As expected, the single pulse widths, 
w s , gleaned from our new objects lie primarily within the 
optimal sensitivity range set by the box-car matched filters 
in our search, as described in section [3] 

The most notable feature of this diagram is that upon 
plotting our integrated pulse width measurements (filled cir- 
cles) with the single pulse width measurements (open cir- 
cles) , it becomes apparent that when using duty cycles mea- 
sured from individual pulses (as were used in the calcula- 
tions of M06, D09, and K09), the values measured for our 
RRAT-like objects become systematically lower by a factor 
of ~3, and up to a factor of more than 10 for the extreme 
case of J0941-39. There may be several causes of this; for 
pulsars with nulling, pulse drifts and mode changes are often 
observed, which can work to broaden the integrated pulsar 
profi le (illustrated well by PSR J2033+00's pulses in Fig. 
|3(d)| above). Multi-component integrated profiles, made up 
of a superposition of otherwise narrow pulses emitted from 
a range in latitudes, can also broaden an object's duty cy- 
cle between single and integrated pulse measurements (as 
has occurred in the extreme case of PSR J0941-39 with 
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Figure 7. The duty cycle versus period for various detections (as 
detailed in £|5.2|l . Among the sources presented in this paper, pul- 
sars are highlighted in blue, RRATs in black, nulling pulsars in 
grey, and PSR J0941— 39 is represented by the red symbols. The 
duty cycles of our detections are shown in two ways: using the in- 
tegrated pulse width (measured after all available on-pulses were 
stacked at the pulsar period) and using the width of the bright- 
est single pulse. M06, D09, and K09 reported only single-pulse 
widths. The dotted lines represent traces of w = 0.125 ms and 
32 ms, beyond which our search sensitivity to individual pulses 
drops significantly. 



10 3 



1(f 



e 

CO 



10" 



10 1 



This work 
D09 
M06 
K09 
„ = 66 mjy 
= 704 mjy 







1 



2 3 4 5 6 7 8 
Distance (kpc) 

Figure 8. Dependence of luminosity on the distance of the 
sources. The lower curve represents our single pulse flux sensitiv- 
ity limit for a pulse width of 32 ms, while the upper curve shows a 
representative limit above which a highly modulated pulsar would 
become detectable in a Fourier search, as detailed in £|5.3I 



its triple-peaked profile). It is clear that were an integrated 
pulse width measured for the sources of M06, D09, and K09, 
the duty cycle distribution of RRATs would more closely re- 
semble that of the general pulsar population, instead of ap- 
pearing to sit characteristically lower tha n pulsars (the com- 
posite profiles for two RRATs shown in iMcLaughlin et all 
(2009) verify this suggestion) . As the RRAT and pulsar duty 
cycle distributions are similar, it is probable that the beam- 
ing fraction of such objects will likewise be akin to that of 
pulsars. 

Next we examine the distribution of the single-pulse lu- 
minosities, as shown in Figure [8] The points cover a broad 



range in both distance and luminosity, appearing to show a 
downturn with proximity to the observer. The lower single- 
pulse sensitivity limit of our search obeying L oc d 2 is plotted 
(at 66 mjy for our 32 ms filter). In light of the possibility 
that scenario (2) (modulation) is true, we consider the fol- 
lowing: Fourier transform sear ches for SNR > 8 pulses as 
done by lEdwards" et all (120011 ) over the 4.4 minute survey 
pointings have a sensitivity limit to pulsars (given a typi- 
cal duty cycle of 5%) with an average single-pulse peak flux 
of (.Speak) *5 3.3 mjy. The prominence and statistics of ex- 
tremely modulated pulsars (like PSR B0656 +14) are not yet 
well ch aracterised, however the study of IWeltevrede et ail 
(2006a|) of 187 pulsars at 1 GHz revealed that typical mod- 
ulation rarely exceeds ten times the average pulse energy 
(with pulses exceeding this range typicall y las ting n anosec- 
onds, and called "giant pulses", ICairn 1 e.g. 120041 ). Their 
follow-up study of PSR B0656+14, as previously noted, re- 
vealed infrequent single pulse emission peaking up to 420 
times the average single pulse intensity. Here we consider 
that such extreme sources will be relatively uncommon, and 
take a population of hypothetical sources bursting at half 
the maximum S ma x/( Speak) ratio of B0656+14, that is, a 
bursting population with typical Smax = 210(S P cak). If such 
a population is contained in our data, sources above the limit 
of (Speak) ^ 3.35 mjy will have bursts at S max — 700 mjy. 
Sources brighter than this limit will be detectable in Fourier 
searches and be tagged as pulsars, however below this limit 
the periodic emission will have been missed, with only the 
brightest bursts detectable in single pulse searches. This 
limit is plotted in Fig. [S] and interestingly, appears to cor- 
respond to the upper envelope of the maximum source lu- 
minosities found by our single pulse search. 

The natural interpretation of this is that the nearby, 
highest-luminosity sources exhibiting erratic emission had 
bright enough low-level emission that instead of being de- 
tected as nulling pulsars or RRATs, they were discovered in 
periodicity searches and taken into the population as pul- 
sars (albeit modulated or nulling ones). When a similar test 
is done for the RRATs found by M06 (whose longer, 35.5 
minute pointings have higher sensitivity to periodic signals 
and a higher likelihood of containing the rarer brightest out- 
bursts of a highly modulated pulsar; a periodicity s earch 
of that survey is described in iManchester et al.l 12001). two 
sources, PSRs J1317-5759 and J1819-1458, still lie signifi- 
cantly above the limit of the Fourier analysis. In the frame- 
work of our interpretation, these two sources are therefore ei- 
ther examples of extreme sources among modulated pulsars 
(however nearby equivalents among the known pulsar popu- 
lation, apart from B0656+14, have yet to be recognised), or 
alternatively they are simply examples of genuine RRATs, 
where the intrinsic pulsar emission is quite strong, however 
usually invisible because of a very high nulling fraction. 

Among the four of our RRATs for which we could mea- 
sure periods, although very statistically poor, the previously 
noted RRAT feature — that they have a tendency to longer 
periods than pulsars — is observed in our new objects, where 
50% have periods longer than one second. Plotted in Figure 
[9] are the normalised distributions for our and the published 
"classic RRAT" sources as compared with the known pulsar 
population for pulsars with P > 50 ms (to restrict the anal- 
ysis to non-recycled pulsars). When a Komolgorov-Smirnov 
test is run on the un-normalised distributions, it results in 
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Figure 9. Normalised period distributions for pulsars and the 
single-burst emitting RRATs of M06, D09, K09, and this survey. 



a high probability, P = 0.999985 (> 4a), that the distri- 
butions are not drawn from the same parent population. 
Few pulsars have periods greater than P — 3 s and because 
sources such as anomalous X-ray pulsars/magnetars/soft 
gamma-ray repeaters all have periods greater than a few 
seconds, it is tempting to draw an association between such 
sources an d RRATs; the recent relea se of timing data on six 
RRATs bv lMcLaughlin et all (|2009l ) show a further similar- 
ity of RRATs to these other neutron star types in their place- 
ment in the period-period derivative diagram, and a simi- 
lar result in periodicity distribution differences. However, 
the broad distribution of source periods and the presence 
of 7 RRATs with periods below one second suggests that a 
more basic factor in the population, for instance age, which 
can create the observed tendency to longer periods. For the 
RRATs which have had timing analysis done, all appear to 
have ages o f a few times 10 5 ~ 6 years. However, the statistical 
analysis of iMcLaughlin et al.l l|2009l ). which compared vari- 
ous derived RRAT and pulsar quantities, does not indicate 
a strong difference in the distribution of RRAT and pulsar 
ages, though does indicate a strong probability that the dis- 
tribution of RRAT and pulsar magnetic field strengths are 
significantly different. There appears to be no relation be- 
tween RRAT period and distance from the galactic plane (a 
proxy, albeit a loose one, for pulsar age). 

The source PSR J0941-39 contributes a peculiar voice 
to this discussion, appearing for hours, to possibly weeks 
at a time like a classic RRAT — erratically emitting single, 
non-sequential pulses at a relatively constant rate — while at 
other times appearing as an intense (~ 10 mjy) pulsar with 
a < 10% nulling fraction, extraordinary sub-pulse drift, and 
latitude-dependent modulation. The object appears super- 
ficially to tie a close link between nulling phenomena and 
RRATs, possibly signifying a pulsar which is in some way 
transitioning from a pulsar mode to an RRAT mode. If this 
is the case, the source provides a link both between RRATs 
and pulsars as well as a link to older, non-radio-emitting 
neutron stars — i.e. that the RRAT phase represents a degra- 
dation of the radio emission mechanism before advancing 
beyond the radio pulsar death line. As a timing solution is 
acquired for this pulsar, it will be interesting to note whether 
it is of an advanced age. It is rather surprising how luminous 
this pulsar is, and what a small nulling fraction it possesses 
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when in its "pulsar" state. It is important to note that the 
site of the RRAT emission is coincident with the three peaks 
of the profile when the pulsar is on. Thus, the RRAT emis- 
sion is from the same site as the normal pulsar emission, 
again implying similar beaming fractions for the two states. 

The properties of the new transients found in this 
search, and particularly of the five new RRATs, give var- 
ious indicators that RRAT behaviour is well described by 
either or both of scenarios (2) or (3) above. The idea that 
observed RRAT behaviour may be an observed consequence 
of extreme (> 99%) nulling is supported by the high nulling 
fraction of the objects tagged as "nulling pulsars" above; 
while previously the "rec ord" nulling fract ion was held by 
PSR B1713-40 at 95% (|Wang et all 120071 ), PSRs J1825- 
33, J1647-36, and J1534-46 all have implied nulling frac- 
tions of more than 97%, telling us that these objects may 
imply a "bridging population" of high-nulling-fraction pul- 
sars, while objects with higher nulling fractions (smaller 
and less frequent on-windows) exhibit RRA T behaviour. A 
full st atisti cal nulling frac t ion an alysis of the lEd wards et all 
(200l|) and Ijacobv et al . (2009) pulsars may indicate that 
the nulling fraction of the RRATs discovered in this search 
fall naturally in the tail o f pulsars with sign ificant nulling 
fractions shown to exist bv lWang et al] (|2007l ) in other sur- 
veys. This proposition, along with the lack of RRAT fea- 
tures that are physically distinct from pulsars (except their 
period distribution) — including their galactic and duty cy- 
cle distribution — point to a strong relationship to the "nor- 
mal" radio pulsar population. The story is complicated 
and strengthened by the apparent mode-switching of PSR 
J0941-39, which lends significant support for scenario (3) 
in that it is clear that the object in its RRAT mode is 
not simply a modulated pulsar, and it is likewise clear that 
the object is both RRAT and pulsar in nature, signifying a 
transition between the two populations and more specula- 
tively, implying that nulling pulsars may undergo stepped 
increases in t heir n ulling fraction with age. The analyses of 
IWang et all (|2007T ) indicate no correlation between nulling 
fraction and characteristic pulsar age, however the number 
of pulsars exhibiting nulls appears to increase with proxim- 
ity to the radio death line (see their Fig. 8). If all RRATs are 
relics of nulling pulsars, the broadening of integrated pulse 
measurements supports the fact that despite their sparse 
emission, the emitted pulses may still be exhibiting drifting 
or mode-changing phenomena. 



6 SUMMARY 

We hav e re-examined t h e inte rmedi ate galactic la t itude sur- 
veys of lEd wards etail i|200ll ) and Ijacobv et ail i|2009h for 
single pulses of radio emission of duration 0.125 < w < 
32 milliseconds. This search revealed 14 intermittent sources 
which appear to be neutron stars. Periods ranging 0.2 to 
6.2 s were measured for 12 of these sources. While only four 
of the objects could be detected retrospectively by Fourier 
search techniques, the remaining 10 exhibited intermittent 
emission behaviour akin to both e xtreme nulling pulsars an d 
the classic RRATs discovered by iMcLaughlin et ail |2006). 
The sensitivity of our search techniques to bursts of lower 
dispersion measure allowed us to find objects which were 
closer than those detected by previous searches. 



The general properties of the new sources appear to 
paint details of the extreme intermittent and RRAT popu- 
lation that concur with having similar emission modes and 
galactic distribution to pulsars, despite being significantly 
less emissive than the general population. It was found that 
the duty cycle distribution of our detections mimics that of 
pulsars. These similarities, and the discovery in this search of 
pulsars with extremely high nulling fractions, lead us to sug- 
gest that RRATs may represent an extreme type of nulling 
pulsar whose window of on-activity is shorter than the rota- 
tion period of the object, and that the sources observed to 
exhibit RRAT behaviour are a natural extension of the tail 
of pulsars observed with high nulling fractions. Furthermore, 
PSR J0941-39 appears to oscillate between a clear "nulling 
pulsar" state and a more intermittent "RRAT" mode. This 
source may represent the first detection of an object in tran- 
sition from a bright pulsar with an average nulling fraction 
to a single-bursting source, closely tying nulling pulsars to 
RRATs and indicating a possible evolutionary progression. 
Further studies and timing are warranted and upcoming on 
this object. 

We conclude by noting that identifying RRATs is not 
a simple task, as they may represent a mixed population 
of intrinsically single-pulsing pulsars and highly modulated 
pulsars whose brightest single pulses singly exceed the sys- 
tem noise. Modulated pulsars appearing as RRATs are ex- 
pected to be few, as among the local population of pulsars, 
PSR B0656+14 is currently the only pulsar known to be as 
modulated as the population required for RRAT emission. 
Nearby pulsars which are part of such a population may in 
the future be identified by pulse energy distribution studies 
of nearby highly modulated pulsars, using long observations 
which have a chance of detecting the rarer brightest out- 
bursts from such objects. Statistical studies of these pulsars 
may help separate "true" RRATs from simply highly mod- 
ulated pulsars. 

To discern the nature of RRATs, other extreme pulsar 
intermittents, and highly modulated objects, it will be crit- 
ical in coming years to both identify more such sources, and 
to continue timing analyses on the known sources to deter- 
mine their relationship to other neutron star populations. 
Of additional interest are the statistics of intermittent phe- 
nomena in extreme nulling/modulated pulsars among the 
known pulsar population, the analysis of which may reveal 
further insight into the relationship between regular pulsars 
and the vast variety of intermittent behaviours observed in 
these equivocal neutron stars. 
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